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Abstract 

Dislocations are at the heart of the plastic behavior of crystalline materials yet it is notoriously 
difficult to perform quantitative, non-intrusive, measurements of their single or collective proper- 
ties. Dislocation density is a critical variable that determines dislocation mobility, strength and 
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ductility. On the one hand, individual dislocations can be probed in detail with transmission elec- 
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tron microscopy. On the other hand, their collective properties must be simulated numerically. 
Here we show that ultrasound technology can be used to measure dislocation density. This devel- 
opment rests on theory — a generalization of the Granato-Liicke theory for the interaction of elastic 
waves with dislocations — and Resonant Ultrasound Spectroscopy (RUS) measurements. The cho- 
sen material is aluminum, to which different dislocation contents were induced through annealing 



and cold rolling processes. The dislocation densities obtained with RUS compare favorably with 
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those inferred from X-ray diffraction, using the modified Williamson-Hall method. 
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I. INTRODUCTION 



First introduced as a curiosity in the mathematical description of elastic continua, disloca- 
tions have become a fundamental building block in the explanation of the plastic behavior of 
crystalline materials. The advent of transmission electron microscopy (TEM) enabled their 
direct visualization and provided a firm experimental foundation for their subsequent study 
Nowadays, dislocations not only describe a physical reality, they also provide a conceptual 
framework to explain very many phenomena in crystal plasticity. However, quantitative 
modeling with predictive power has remained elusive. 

At the atomistic level, there has been significant progress in the theoretical, direct obser- 
vation, and numerical modeling of individual dislocations. A few examples are the work of 
Leyton et al. [1] on the prediction of energy barriers to dislocation motion, the direct obser- 
vation of dislocation mechanisms of Oh et al. [2], and the molecular dynamics simulations of 
Marian et al. [3]. However, in spite of the considerable advances in the study of individual 
dislocations, the fact remains that the collective behavior of many dislocations, a critical 
ingredient in crystal plasticity, is not necessarily obtainable by simple aggregation of the 
individual behavior of many isolated dislocations. In fatigue, for example, a non-intrusive 
characterization of the material prior to fatigue failure would represent a significant step 
towards the control of fatigue damage. Progress in this respect has been achieved through 
the use of infrared thermography [3HS] - On the numerical modelling front, a direct approach 
to tackle the problem of the collective behavior of dislocations is the Dislocation Dynamics 
Method [7HTU]. In a related development, Warner et al. [H] have performed a multiscale 
simulation that sheds light on the influence of dislocation formation on the propagation of 
cracks in aluminum. 

Dislocation density is closely related to the ability of a given material to endure plastic 
deformation [T2HII]- But, can dislocation density be accurately measured? TEM allows for 
a local measurement in very small specimens that have to be especially prepared, including a 
non direct thickness measurement to estimate the local volume of the area under observation 
[15] . This technique is of limited use when individual dislocations cannot be resolved due 
to high densities |16j . Another widely used possibility is to consider the broadening of X- 
ray diffraction (XRD) peaks by dislocations and other defects [IT], a formulation that has 
been applied to submicron grain size and heavily deformed copper [TTH2T)] . nanocrystalline 
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iron powder and iron-based powders [2TH23] . and aluminum alloys [24H26] . Recently, Balogh, 
Capolungo and Tome [TB] have provided an assessment of this method, including its strengths 
and limitations. Clearly, an alternative, reliable, non-intrusive, way to measure dislocation 
density would appear to be a desirable development. 

Ultrasound (US) technology has been widely used as a non-destructive evaluation tool for 
several decades |27J. This is based on the fact that the associated ultrasonic energy is very 
low compared to energies needed to deform, break, or significantly alter the tested material 
in any way and on the further fact that interfaces, as well as flaws, within a material, affect 
the propagation of acoustic waves. So, monitoring ultrasound transmission or reflection 
provides information about the insides of a material or structure. 

Dislocations interact with elastic waves. This has been known since the early days 
of dislocation theory [28| 129] . Can ultrasound technology be used to learn something — 
quantitatively but non-intrusively — about dislocations in a crystalline solid? The effect that 
dislocations have on nonlinear wave propagation has been studied [30, 31J, as well as the 
possibility of using it to characterize fatigue microstructures [32] • More specifically, in re- 
cent years Resonant Ultrasound Spectroscopy (RUS) has emerged as a very efficient tool 
for materials characterization [3"3Tl3~T] . A complete set of elastic constants can be measured 
precisely for any kind of symmetry. It has proved to be useful in different fields, as in 
condensed matter where its precision has been essential in order to demonstrate that a 
quasi-periodic lattice is indeed elastically isotropic [3S]. Other examples are the observation 
of giant softening in high temperature superconductors [39J, the measurement of the me- 
chanical properties of biomaterials [40], the determination of elastic constants of thin films 
[41] . the mechanical characterization of composite materials [42], the measurement of inter- 
nal friction in polycrystalline copper [43], the characterization of elastic anisotropy in a bulk 
metallic glass [44], and its application in geophysics, in particular for the measurement of 
elastic constants at high temperature [43] . A recent publication [3H] compares RUS to other 
more traditional techniques, including four-point bending, nanoindentation and impulse ex- 
citation, concluding that it is much more precise than most of them for the measurement of 
elastic constants. 

In this article we show that, other things being equal, RUS can distinguish between 
materials with different dislocation densities. This is demonstrated preparing samples with 



different dislocation content from the same as-received aluminum bar (Section II A), and 
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measuring their elastic constants and elastic wave velocity with RUS (Section II B). This 
provides a measure of their dislocation density (Section |lII[ ). The results are further compared 
and corroborated with those obtained by XRD peak broadening profile analysis (Section IV). 



II. EXPERIMENTAL 
A. Sample Preparation 

Commercially 1100 pure aluminum (99.0% pure) was used to perform RUS and XRD 
measurements. From the same as-received bar, five pieces were taken to prepare the studied 
conditions classified as original, annealed and laminated material: two samples were annealed 
at 673 K, one for 5 h and another for 10 h. Two others were cold-rolled, at either 33% 
or 43% of the initial diameter of the as-received original material. It is well known that 
longer annealing leads to lower dislocation density, and stronger cold-rolling leads to higher 
dislocation density. However, there is no known way to quantitatively estimate how much 
the dislocation density will change as a consequence of a given amount of annealing or cold- 
rolling. The pieces were thus numbered from 1 to 5, in the sense of expected increasing 
dislocation density. From each one of the five pieces, presumably with different dislocation 
densities, one portion was set aside for RUS testing, and another for XRD. 

The five RUS samples were shaped as rectangular parallelepipeds, with dimensions pre- 
sented in Table [TJ Opposite sides are parallel within 0.06° and adjacent sides are orthogonal 
within 0.3°. Our samples can be then modeled as perfect rectangular parallelepipeds [l6"l W7\ . 
Samples are labeled 1,2,3,4 and 5 for increasing expected dislocation density. 

All conditions tend to have anisotropic grain size distribution due to the conformation 
by extrusion in the case of original (as-received) and annealed bars, and due to lamination 
in the case of laminated conditions. Considering this, all samples for XRD were cut in the 
transversal axial direction from the original and annealed bars, and transversal to the lam- 
inated direction in laminated conditions. XRD samples were further chemically attacked 
with a solution of 30% HC1, 10% HF and 60% distilled water to remove superficial defor- 
mation induced during sample cut. A Siemens D5000 diffractometer was used for XRD 
measurements, using Cu Ka (1.5418 A)radiation. The instrument-broadening contributions 
were measured independently using (111) and (100) oriented silicon single crystals, and were 
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Sample 1 


Sample 2 


Sample 3 


Sample 4 


Sample 5 




Annealed 


Annealed 


Original 


Rolled 


Rolled 




673 K /10 h 


673 K/5 h 




at 33% 


at 43% 


d\ (cm) 

J- \ / 


1.701 ±0.001 


1.700 ±0.001 


1.7011 ±0.0002 


1.696 ±0.001 


1.7007 ±0.0006 


d 2 (cm) 


1.0015 ±0.0002 


0.9997 ±0.0004 


0.9998 ± 0.0003 


1.001 ±0.001 


1.001 ±0.001 


d 3 (cm) 


4.902 ±0.001 


4.900 ±0.001 


4.901 ± 0.001 


4.901 ±0.001 


4.900 ±0.001 


M (g) 


22.45 ±0.01 


22.38 ±0.01 


22.43 ±0.01 


22.35 ±0.01 


22.42 ±0.01 


p (g/cm 3 ) 


2.688 ±0.002 


2.687 ±0.002 


2.691 ±0.002 


2.687 ±0.004 


2.687 ±0.003 



TABLE I: Sample characteristics: Rectangular parallelepiped dimensions, mass and mass den- 
sity for the five samples. Columns are ordered for increasing expected dislocation density. 

subtracted from the data prior to calculations [18] . Si(lll) is located at an angle of 28.426° 
while Si(400) is at 69.192°, with a broadening measure at the Full Width at Half Maximum 
(FWHM) of 0.022° and 0.034°, respectively. Si(lll) FWHM was subtracted in the case 
of the Al(lll) and Al(200) signals while Si(400) was used for Al(311) and Al(222). Three 
XRD measurements were performed for samples 1, 2, 4 and 5, while two measurements were 
performed for sample 3. 

B. Resonant Ultrasound Spectroscopy 

An in- house built RUS apparatus |I9J [50] , similar to those presented in [331 G2] was used to 
measure the speed of elastic waves of the five aluminum samples. This was done through the 
measurement of the resonant frequencies of the samples for stress-free boundary conditions. 
These frequencies, in turn, provide the elastic constants (in Voigt notation) after the 
mass and linear dimensions of the samples have been independently measured. In order to 
satisfy this boundary condition, samples are placed in between an emitter and a receiver, 
held by two opposite corners, the receiver being held by a set of springs mounted on a linear 
air bearing. The weight of the receiver part of the setup is such that, at equilibrium with the 
spring force, the distance between the emitter and receiver surfaces is slightly larger than 
the sample's diagonal. In this way, the sample can only be held by its corners in the setup 
if a small mass is added to the receiver part, typically of 10 g. Thus, the sample- apparatus 
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contact force is small, with an upper bound of 0.1 N, corresponding to about one half of the 
sample weight Mg ~ 0.22 N. The estimated contact area is ~ 0.1 — 0.2 mm 2 . 

The results of RUS measurements are shown in Table [II] and Figure[TJ Preliminary results, 
with a smaller data base, were reported in [51]. A satisfactory fit to the data is obtained 
under the assumption of transverse isotropy, so that, of the 21 possible elastic constants, 
only five are independent. C44 can be determined with much higher accuracy than the other 
elastic constants, a well known characteristic of RUS [37]. From it, and the independently 
measured mass densities, the speed of shear waves for each sample can be determined with 
an accuracy of ~ 0.1%. 

For each one of our five samples a set of 30 resonant frequencies was measured, starting 
from the lowest resonant mode around 19 kHz and finishing at ~ 150 kHz. The pure elastic, 
non-dissipative assumption is well verified: the lowest quality factor Q is ~ 10 3 . In order to 
verify that all resonance frequencies are measured in the elastic linear regime, each resonant 
curve is measured for five different driving amplitudes. We then verify that the resonance 
amplitude is linear with the driving amplitude. We emphasize that the main difference with 
our previous, preliminary, results presented in [51] is the number of independent measure- 
ments for each sample. For the results presented in the present manuscript each sample 
was placed ten times in the apparatus in order to reduce errors due to slight dependence of 
the resonant frequencies on the contact load and positioning with respect to the ultrasonic 
receiver. In [51], each sample spectrum was measured just once. In fact, there is an incorrect 
phrase in [5T] that states that each sample was measured five times. This was the case in 
early measurements for a lower number of measured frequencies (12) assuming complete 
isotropy. 

As briefly reported in [51], attempts to fit resonant frequencies assuming homogeneity 
and isotropy fail considerably. In fact, under the assumption of isotropy, measurements of 
C44 agree with previous published results but not those of Cn, which are about 20% lower 
than the values reported in the literature. This discrepancy arises independent of sample 
preparation, positioning and the particular emitting transducer type. The complete isotropy 
hypothesis is then questionable. Indeed, our aluminum samples are prepared with their 
longest dimension in the extrusion direction of the original aluminum bar and microscopic 
images (not shown) confirm that crystal grains are elongated in this direction. We then 
consider that our samples are transversely isotropic, which implies that the elastic-constant 
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Sample 1 


Sample 2 


Sample 3 


Sample 4 


Sample 5 






Annealed 


Annealed 


Original 


Rolled 


Rolled 






673 K /10 h 


673 K /5 h 




at 33% 


at 43% 


RUS 


C u (GPa) 


111 ± 1 


112 ± 1 


110.2 ±0.3 


112 ± 1 


112 ± 1 




C33 (GPa) 


112 ± 1 


113 ± 1 


111 ± 1 


113 ±3 


112 ± 1 




C23 (GPa) 


61 ± 1 


62 ±1 


60.5 ±0.4 


63 ±2 


61 ± 1 




C12 (GPa) 


58 ± 1 


60 ±1 


58.2 ±0.3 


60 ±2 


60 ±1 




C44 (GPa) 

'±'± \ / 


27.52 ±0.08 


27.34 ±0.03 


27.40 ±0.02 


27.11 ±0.06 


27.02 ± 0.03 




v L (m s" 1 ) 


6412 ± 19 


6467 ± 22 


6399 ± 9 


6450 ± 35 


6450 ± 21 




t"r (m s" 1 ) 


3200 ± 5 


3190 ± 2 


3191 ± 1 


3177 ±4 


3171 ±2 


XRD 


A(10 8 mm -2 ) 


11.4 ±0.5 


9.1 ±4.5 


14.3 ±8.9 


63 ±8 


65 ±6 




D (nm) 


86 ±27 


82 ± 17 


94 ±22 


130 ± 34 


120 ± 36 



TABLE II: RUS and XRD measurements results. Top seven lines: Elastic constants (Cu, C33, 
C23, C12 and C44 in Voigt notation), longitudinal wave velocity (vl), and shear wave velocity (vt) 
obtained through RUS. Columns are ordered for increasing expected dislocation density. Absolute 
errors for Cij constants are computed from the standard deviations of ten RUS measurements 
obtained for each sample. The errors of vl and vt are obtained from the standard deviations 
of the set of ten values obtained from their definitions vl = y/Cu/p and vt = \fC~u]~p, with p 
the mass density. For transverse isotropy, the other elastic constants are C22 = Cu, C13 = C23, 
C55 = C44 and 2C66 = Cu — C\i- Note that C44 can be determined with much better accuracy 
than the other constants. Bottom two lines: Dislocation density A and crystallite size D, obtained 
by XRD. Absolute errors for A and D are computed from the standard deviations of three XRD 
measurements performed for samples 1, 2, 3 and 4, and for two measurements performed for sample 
3. 

matrix has five independent values. Having that in consideration the measured resonant 
frequencies fit much better the theoretical predictions. It is important to stress that this 
is not just due to the increment in the number of adjustable parameters, it is truly a 
consequence of a better assumption concerning the material structure. To prove this, we have 
performed anisotropy tests by varying the anisotropic and isotropic planes in the RUS inverse 
code. The better results, quantified by the final root-mean-square error between measured 
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3195 
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3180 
3175 
3170 
3165 
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2 3 4 

Sample number 



FIG. 1: Elastic constants and wave velocities of samples with different dislocation 
content as measured by RUS. Elastic constants C\\ and C44 (top) and longitudinal {v£) and 
transverse (vt) (bottom) sound speeds as determined by RUS for the different aluminum samples. 
Open circles correspond to results obtained for ten different sample positionings and solid circles 
correspond to the mean values. Error bars are computed from standard deviations of the set of ten 
measurements. Note the difference in scale between the left-hand-side-panels and the right-hand- 
side ones. The shear modulus, C44, is determined with much higher accuracy than Cn, allowing 
the speed of shear waves, vt, to be determined with an accuracy of ~ 0.1%. It exhibits a clear 
downward trend in the direction of expected increasing dislocation density. 

and predicted frequencies, correspond to those where the sample's longest dimension is 
the anisotropic one with respect to the transverse isotropic dimensions. The measured 
elastic constants, as well as the longitudinal (vl) and shear (vt) wave velocities related to 
Cn and C 44 respectively, are given in Table I. Figure 1 shows Cn, C 44 , Vl and versus 
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sample number. Each elastic constant is obtained by averaging 10 results obtained from 
10 independent sample positionings. The reported errors correspond to standard deviations 
of these 10 measurements. As expected, deviations from isotropy are small, for example 
C11/C33 > 0.99 and C^/Cqq < 1.06. Thus, a comparison with the theory presented below is 
possible. No clear tendency can be observed for vl because of the large experimental errors, 
i.e., 15%. However, vt data is much more precise, errors ~ 0.03% — 0.15%, and a clear 
decreasing tendency in the direction of expected increasing dislocation density is observed. 

III. INTERPRETATION OF RUS MEASUREMENTS 

In order to link the RUS measurement results to dislocation density we use the formula 



Eqn. ([I]) links the relative change in shear wave velocity Avt/vt between two samples of a 
material that differ in dislocation density n, where n is the number of dislocation segments 
of length L, Burgers vector b and line tension T, per unit volume. The subscript e is for 
edge dislocations, s is for screws, and \x is the shear modulus of the reference material. This 
is a result valid to leading order in perturbation theory when the changes in vt are small. 
Formula (ITJ) shows that an increasing dislocation density results in a decreasing speed of 
shear waves. It is obtained using Multiple Scattering Theory, a reasoning that quantifies 
the intuition that dislocations will make it more difficult for waves to propagate [31]. The 
physical process that is responsible for the change in the speed of shear waves is as follows: 
an elastic wave is incident upon a dislocation segment of length L whose ends are pinned. 
As a result, it oscillates like a vibrating string. This oscillation, in turn, generates secondary 
waves. As explained in Appendix [A] the coherent superposition of many such scattering 
processes leads to Eqn. (II]). 
Considering that [33J 



[321 E3] (see Appendix [A} 




(1) 




(2) 



(3) 
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where vl is the speed of longitudinal waves, we have 

^ « ~< + = -^jA(nL") . (4) 

What RUS does, then, is to provide a measure of the dimensionless quantity nL 3 . This 
formula does not provide an absolute measurement of dislocation density, but a measurement 
of the difference in dislocation density between two samples. The data show that our RUS 
measurements cannot, within experimental error, differentiate between samples 2 and 3, but 
can differentiate sample 1 from sample 2, sample 3 from sample 4, and sample 4 from sample 



5, providing the values of dislocation densities given in Table III We have (Avt/vt) ~ 
A(nL 3 )/60 with the values given in Table III Of course, the same value of A(nL 3 ) can be 



obtained with a variety of values for n (number of segments per unit volume) and L (distance 
between pinning points) separately Taking, for visualization purposes L ~ 10 nm, we get 
that the samples differ in dislocation density AA ~ (1 — 3) x 10 9 mm -2 (in the usual units). 





Samples 1-2 


Samples 2-3 


Samples 3-4 


Samples 4-5 


x 10 3 

(v T ) 

A(nL 3 ) 


3.1± 1.7 
0.19 ± 0.11 


- 0.3 ± 0.7 
- 0.02 ± 0.04 


4.4 ± 1.3 
0.28 ± 0.08 


1.9 ± 1.4 
0.12 ± 0.09 



TABLE III: RUS measurement of dislocation density. RUS measures the difference Avt in 
the speed of shear waves between the various samples. Results are given in the first line relative 
to the arithmetic mean (vt) of each pair of values. From these values and using Eqn. Q, the 
difference in dislocation density A(nL 3 ) can be obtained (second line). Within experimental error, 
it is not possible to distinguish the difference in dislocation density between samples 2 and 3. 

We now provide validation to the dislocation density values obtained through RUS with 
XRD peak-broadening measurements on the same samples. 



IV. X-RAY DIFFRACTION ANALYSIS 



A typical XRD pattern obtained from the aluminum samples is shown in Figure 2. The 
characteristic (111), (200), (311) and (222) signals are observed while the (220) reflection is 
absent. The absence of the (220) reflection is due to texturing induced by the deformation 
of the sample, a well known phenomenon in metallurgy |56j. For all peaks present in the 
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diffractograms we observed that the larger the matrix deformation (from the annealed to 
the extreme laminated conditions) the larger the peak width. The inset in figure 2 shows 
this effect for the Al(lll) reflection. 




40 50 60 70 80 

29 (°) 

FIG. 2: XRD pattern for Sample 2. The characteristic (111), (200), (311) and (222) signals 
are observed while the (220) reflection is absent due to texturing induced by the deformation of 
the sample. The inset shows a zoom of the (111) peak of the XRD patterns for the five samples. 

Current theory [57] accounts for the broadening of XRD peaks as arising from two effects, 
finite crystallite size and presence of dislocations. If the two effects are uncorrelated, they 
will act additively: 

AK = AK S + AK D (5) 

where AK is the Full Width at Half Maximum (FWHM) of the diffraction peak at wave 
vector K. AK is the contribution of the crystallites finite size, and AK D is the contribution 
of the dislocations. On dimensional grounds AK S — C\jD where D (dimensions length) is 
the crystallite size, and C\ is a dimensionless constant. This quantity is independent of the 
X-ray wavelength. 

Let us estimate AK D on dimensional grounds [58]. It must depend on b = \b\, the mag- 
nitude of the Burgers vector of the dislocations (assuming only one value for its magnitude 
and taking an average if there are several; in any case its value is of the order of the size 
of the cube root of the volume of the unit cell). Also, it must depend on A, the dislocation 
density (dimensions length -2 ), assumed uniform. Consequently, it must involve the prod- 
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uct 6A 1 / 2 , since the broadening must be an increasing function of b (larger b means greater 
deformation, hence larger broadening). For low dislocation densities (i.e. such that AK 
is small compared to K), its dependence must be linear. There is an additional quantity 
with dimensions of inverse length, K, the diffraction vector of the broadened peak, and for 
dimensional analysis consistency we get 

AK D = C^ /2 bA 1/2 K . (6) 

There cannot be an added contribution from the dimensionless quantity bK since in that 
case there would be a broadening even for vanishing dislocation density A. Here C2 is a 
dimensionless constant that may depend on the crystalline plane (Miller indices (h,k,l)) 
doing the scattering, but it does not depend on A. Higher powers of b (with concomitant 
higher powers of K) seem unlikely, again claiming linearity at the lowest order. This type 
of functionality is deduced through the modified Williamson-Hall method (59j [60] from an 
analysis of the full shape of the spectral lines (see Appendix [B]) , with the result that 

where C^oo an d q can be expressed in terms of the elastic constants of the crystal and 
empirically determined parameters, H 2 is the fourth order invariant of the hkl indices of the 
different reflections and M 2 is a constant depending on the effective outer cut-off radius of 
dislocations. We take M 2 = 0.1 [59j EQ]. Consequently, a plot of (Ai^) 2 vs K 2 will yield 
a straight line: The slope determines the dislocation density A, and the intercept with the 
vertical determines the crystallite size D. Eqn. Q deals with (AK) 2 , the square of the 
FWHM. The result is equivalent to the result for AK obtained in ^ using dimensional 
analysis as long as crystallite size and dislocation density are uncorrelated [2~T] . 



A. Modified Williamson-Hall plot 

Figure [3] shows the (AK) 2 vs. K 2 C plot calculated according to the 'modified' 
Williamson-Hall method from the values obtained of the diffraction patterns following Eqn. 
0. The increasing slope from annealed to laminated conditions indicates an increasing 
dislocation density according to Eqn. ([TJ). The values of dislocation density A for all condi- 
tions are shown in Table [XTJ The intercepts with the vertical axis are, within experimental 
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error, similar, indicating no significant change in the crystallite size D among the different 
samples, so that the X-ray peak broadening is due only to the dislocations (Table |H|. 
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FIG. 3: AK 2 versus K 2 C plot for the five different aluminum samples, obtained through 
the modified Williamson-Hall method. According to Eqn. ([7]) the slope provides a measure of 
the dislocation density, and the intercept with the vertical axis provides a measure of the crystallite 
size. The resulting values are quoted in Table [TIJ 



V. CONCLUSION 

Qualitatively, Figure [T] shows that, as expected, the shear modulus C44 of the different 
samples decreases in the direction of expected dislocation density. Also, Figure [3] shows 
that XRD lines also broaden, as expected, in the direction of expected dislocation density. 
Figures [T] and [3] show that the five samples can be separated in two groups: Samples 1, 
2 and 3 with low dislocation density, and samples 4 and 5 with clearly large dislocation 
density. Interestingly, ultrasound and XRD measurements coincide in their diagnostic of 
differentiating the two groups and which samples have the larger number of dislocations. 

Figure [|provides a comparison between the results obtained with RUS and those obtained 
with XRD. Quantitatively, RUS is a better discriminant of dislocation density among these 
five aluminum samples than XRD: The latter can distinguish, within experimental error, 
two different values, while the former can distinguish four. Now, "dislocation density" has a 
slightly different meaning in the two contexts: In XRD, the measurements are modelled in 
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5 6 

A (x 10 9 mm -2 ) 



FIG. 4: RUS and XRD characterization, compared and contrasted. Shear wave velocity 
vt determined through RUS, plotted against the dislocation density A determined through XRD, 
for the five aluminum samples. RUS coincides with XRD in identifying two groups: samples 1, 2 
3, and samples 4, 5. However, RUS goes beyond XRD since it can differentiate 1 from 2 and 3, 2 
from 4, and 4 from 5. The resulting values for A(nL 3 ) are quoted in Table 
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terms of A, with dimensions inverse length squared, denned as the total length of dislocation 
line within a sample divided by sample volume. In RUS, measurements are modelled in terms 
of nL 3 , with L the length of the vibrating dislocation line segments, or distance between 
dislocation pinning points, and n the number of such segments per unit volume. Although 
the theory of [53] is sufficiently general to accomodate dislocation segments of different 
length, here we take all segments of the same, average, length for simplicity. The two 
values of A detected by XRD are ~ 10 9 mm -2 (samples 1, 2 and 3) and ~ 6 x 10 9 mm -2 
(samples 4 and 5). RUS does not provide absolute values of nL 3 , but differences between 
different samples. The precise numbers depend on the nature (screw, edge, prismatic) of 
the dislocations. In our case we can discriminate values of A(nL 3 ) in the range 0.1 — 0.3. If 
we take, tentatively, A = nL and the values of A provided by XRD, we can conclude that 
samples 1 and 2 have the same A but different L, and the same for the samples 4 and 5. 
But samples 3 and 4 differ both in n and L. 

Finally, our XRD measurements show a texture effect due to the different samples history, 
that could lead to missing part of the dislocation density information. Eventually, this could 
be avoided taking cuts in the longitudinal direction of the material for the different conditions 
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[25] . In this sense, RUS presents an advantage over XRD because it takes an average over 
the whole sample in every measurement. 

In conclusion, our results clearly establish that RUS can distinguish between aluminum 
samples with different dislocation densities, thus opening the way for the use of ultrasonic 
techniques in the characterization of the plastic behavior of materials. 
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Appendix A: Interaction of elastic waves with dislocations 

The theory that we use to interpret the RUS measurements is based on recent work by 
Maurel et al. in which the scattering of an elastic wave of arbitrary polarization by a pinned 
(edge) dislocation segment \§T\ of finite length L, or an individual prismatic dislocation 
loop [02] was calculated in detail. The ensueing results were used within the framework of 
multiple scattering theory to develop formulae for the effective wave velocity, and attenu- 
ation, of an elastic wave propagating within an elastic medium filled with edge dislocation 
segments, or prismatic dislocation loops [62J. The case of screw dislocations is calculated 
below. These developments build on earlier results of Granato and Lucke [631 El] and their 
primary advantage for the purpose of the present work is the ability, absent in [B"3"l IB"!] , to 
distinguish between longitudinal and transverse waves. This ability is essential to interpret 
the RUS results. 

Consider an homogeneous, isotropic, elastic continuum, of density p, shear modulus fi 
and bulk modulus B. This is the reference material. It does not matter whether it does or 
does not have dislocations to begin with. All that matters is that it can be described as ho- 
mogeneous, continuous and elastic. Next, take the same medium with an additional random 
(uniform) distribution of dislocations: density (number per unit volume) n e (resp. n s ) of 
pinned edge (resp. screw) dislocation segments of length L e (resp. L s ) and Burgers vector 
b e (resp. b s ). If the material is loaded by a monochromatic elastic wave, the dislocations 
will interact with the wave and will influence the material's response. The model used in 
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all cases is the string model of Koehler [65] : The dislocation is treated as an infinitely thin 
string endowed with mass and elastic tension; as it is loaded by the elastic wave it responds 
by oscillating along glide planes and, as it oscillates, it re-radiates. This process leads to a 
new ( "renormalized" ) value fi R for the shear modulus. The bulk modulus does not change. 
The result is, in the limit of low frequencies, 

J- 1 --^ t F n ' L --^ t F n ' L " (A1) 

from which Eqn. (jl]) follows. An equation, qualitatively similar to ( Al ) that links dislocation 



density with shear modulus decrement, has been provided in [30J. The formulae of [30] 
do provide the correct scaling of said decrement in terms of dislocation length. However, 
they depend on undetermined multiplicative parameters, the resolving shear factor and the 
conversion factor from shear strain to longitudinal strain. The reasoning below provides a 



justification of the numerical coefficients presented in (Al). It also provides a justification 



of the result of [30j as the low-frequency limit of a finite-frequency formulation. 



We shall ignore the effect of prismatic dislocation loops. We now derive the result (Al) 
for screw dislocations. The reasoning follows closely an analogous calculation of [51] for edge 
dislocations. 

Elastic waves in an homogeneous, isotropic, three-dimensional, infinite elastic medium 
are described by displacements u(x, t) as a function of equilibrium position x and time t 
that satisfy the wave equation 

d 2 Ui d 2 u k 
P ^- C ^^ = ° (A2) 

with Cijki = X5ij5ki + fi(5ikSji + 5uSjk) the tensor of elastic constants, and i,j,k,= 1,2,3 
which in the isotropic case has only two independent constants: /i, and A = [B — (2/x)/3]. 
Such a solid supports two types of propagating waves: longitudinal (acoustic) and transverse 
(shear) waves with propagation velocity vl = ^(A + 2/i)/p and vt = yjijp, respectively. 
Their ratio 7 = Vl/vt is always greater than one, i.e. 7 > 1. 

Dislocations are modeled as "strings" of length L, mathematically described through a 
position vector X(s,t), with s a Lagrangean parameter to label points along the string, 
which is pinned at the ends. Their simplest equilibrium position is a straight line. They are 
characterized by a Burgers vector b, perpendicular to the equilibrium line for an edge and 
parallel to it for a screw. Their unforced motion is described by a conventional vibrating 
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string equation 

d 2 X dX t d 2 X t 

m ^ + B ^r- T ^ = (A3) 

where the mass per unit length m and line tension V for screw dislocations are given by 

- -^InfA). (A4) 



4tt V.V 

5 and 5q are external and internal cut-offs, and B is phenomeno logical drag coefficient. 
We shall only consider glide motion since climb involves diffusion, which happens, at room 
temperature, on a time scale long compared to the inverse of the wave frequencies under 
consideration. 



Eqns. ( A2 ) and (A3 ) describe waves and dislocations that do not interact. The interaction 
is introduced through right hand side — source — terms. Also, it is convenient to describe the 
waves not through particle displacement u but in terms of particle velocity v = du/dt. When 
N dislocations are present this leads to 

d\ d 2 v k 

p m^- c - kl d^r Sl (A5) 

where the right-hand-side term Sj is given by 

JV 



Si(x,t) = Cjjkitmnk / ds X^(s,t)r n bi 



n=l 

x^-5(x-X n (s,t)). (A6) 
dxj 

Here, e mn k is the completely antisymmetric tensor of order three, f is a unit tangent along 
the dislocation line. For screw dislocations in equilibrium, the Burgers vector points along 



this tangent: b = bf. In the case of the string equation (A3) the coupling with elastic waves 
is provided by the Peach-Koehler fBS] force : 



mX k (s, t) + BX k (s, t) - rX' k '(s, t) = fib N kjp d 3 u p (X, t), (A7) 

with N k j p = e k j m T m r p + e kpm r m Tj . Overdots mean time derivatives, and primes mean deriva- 
tives with respect to s. 



The procedure now is as follows: The loaded string equation (A7) is solved in terms of 



normal modes, and the solution plugged into the right hand side of the wave equation (A5) 
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In the long wavelength limit, A ^> L (for L up to 100 nm this allows for frequencies well 
into the hundreds of MHz regime), and for small string displacements, the result of this 
operation is 

V ik v k (A8) 



2 d v k 

PU Vi - Cijkl 



where 



N 



dxjdxi 




(A9) 

7r* m W ' " ... 

with S(oj) ~ u 2 /(u 2 — u\ + iuB/m) and u\ is the frequency of the fundamental mode of 
the string with fixed ends: uj\ = (jr/L)y/T/m. 

The next step is to consider wavelengths long compared to the mean distance between 
dislocations, an assumption that is well satisfied by our experimental conditions involving 



centimeter-sized samples, so that the discrete sum in (A9) is smoothed, replacing it with 
an integral over space with a continuous density n(x) of dislocation segments, and the 
tensor N^N^. by its angular average, (N^N^ fe ), assuming all directions equally likely. A 
reasoning like this one is used to study waves in plasmas [67J and it is valid for wavelengths 
long compared to inter-dislocation distance. It is straightforward to check that 



(A10) 



Eqn. (A8) thus becomes, in the case of uniform dislocation density n(x) 



n, 



2 d 2 v k 
-fXJ Vi - c ijk r 



A 



' dxjdxi 
4 2 

^ij8 k i + - {5 ik 5 j i + 5u5 jk ) 



d 2 v k 
dxjdxi 



(All) 



where 



A = (bfj)' 



8 L 



it 2 m uj 2 



-n 



+ iuB/m 

From this expression it is easy to read the renormalized values of the Lame constants: 



A 



R 



R 



2 . 



(A12) 



We see that 3 A + 2p R = 3A + 2/x. Finally, for low frequencies, so u u% and u B/m, 



we get the first term on the right-hand-side of Eqn. (Al). 
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Appendix B: XRD theory: Peak broadening due to the presence of dislocations 



Using quantitative models it is possible to correlate microstructural parameters such as 
crystallite size and dislocation density, with the broadening of XRD peaks induced by stress 
of the crystalline structure during deformation [6EJ EH] . Williamson and Hall [70] suggested 
that the two contributions to the broadening can be represented as 



where AK = 2 cos 6 (AO) /A is the full width at half maximum (in radians), 6 is the diffraction 
angle, A is the wavelength, D is the crystallite size and AK D is the strain contribution to 
line broadening. The Williamson-Hall plot of AK as a function of the scattering vector 
K = 2 sm9/X was modified by Ungar and Borbely [59], replacing K or K 2 by KC 1 ' 2 orK 2 C, 
respectively. These latter procedures, called 'modified' Williamson-Hall method, can be used 
in the case of systems of isotropic particle sizes for a physically correct determination of the 
apparent size and mean square strain [TTHT5] . Thus, if dislocations are the main contributors 
to the residual strain, the modified Williamson-Hall plot can be expressed as: 



where M is a constant depending on the effective outer cut-off radius of dislocations, and 
H 2 = (h 2 k 2 + h 2 + l 2 + k 2 l 2 )/(h 2 + k 2 + I 2 ) 2 for a cubic crystal system. C h00 is the contrast 
factor for Bragg reflection (hOO), which determines the contrast factor for Bragg reflection 
(hkl) introduced by Ungar and co-workers to take into account the influence of the residual 
strain over the different Bragg reflections [211 EOj E], through C = CV)o(l ~ iH 2 ). From 
the linear regression of the left-hand side of Eqn. ^ versus if 2 , the parameter q can be 
determined experimentally. Subsequently, the dislocation density A can also be determined 
from XRD results. 
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